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STABILITY OF THE LAMINAR BOUNDARY LAYER IN A STREAMWISE CORNER 


By 

Wi 1 1 i am D. Lakin* 

SUMMARY 

Viscous incompressible flow along a streamwise rectangular corner form- 
ed by the intersection of two seni -inf inite flat plates is often called the 
corner boundary layer problem. Theoretical work on the mean flow in this 
problem has primarily been asymptotic and numerical in nature. Unfortunate- 
ly, experimental results do not recover the flow patterns predicted by 
theory and, indeed, tend to differ from study to study, Zamir (1981) has 
suggested that one reason for the differences is an instability of the mean 
flow. The present work is the first theoretical investigation of the sta- 
bility of the corner boundary layer problem. 

By symmetry, it is sufficient to consider the flow region y < z in 
which X, y, and z are respectively the streamwise, vertical, and spaawise 
coordinates, and the plate is at y = 0. Rubin (1966) has shown that for x 
large and fixed and z/x small, the presence of the corner induces a secon- 
dary spanwise velocity w which reverses direction within the boundary lay- 
er. Further, it can be shown that with z/x small, there is a wide range 
of z such that, to suitable order, all three velocity components u, v, 
and w are functions of y alone, so the mean flow is quasi -paral lei . 

This region will be termed the blending boundary layer. 

In the present work, three-dimensional, time dependent, small amplitude 

*Eminent TrWel^?oy,' ‘Department of Mathematical Sciences, Old Dominion 
University, Norfolk, Virginia 23508. 
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perturbations u , v , w , and p are superimposed on the mean velocity 
and pressure fields in the blending boundary layer region. As weak stream- 
wise vorticity perturbations may produce large defects in the mean stream- 
wise flow, two different velocity scales are used for the velocity perturb- 
ations. This is similar to the treabnent of Gortler vortices by Floryan and 
Saric (1982). Partial differential equations for disturbance quantities are 
obtained by substituting total velocity and pressure fields into the Navier- 
Stokes equations and consistently linearizing about the mean flow. In ac- 
cord with the quasi-parallel nature of the mean flow, dependence on x, z, 
and t is separated out by assuming normal mode forms for the perturba- 
tions. In the Gortler context. Hall (1983) has shown that there is a wave 
number regime for which parallel flow gives a self-consistent approximation 
to the linear stability equations. 

After soma manipulation, the coupled ordinary differential equations 
for the perturbation amplitude functions i3(y), v(y), w(y), and p(y) can 
be reduced to a main uncoupled differential equation for v, a second dif- 
ferential equation which couples w to v, and algebraic equations for u 
and p. Some care must be taken in formulating associated 'outer' boundary 
conditions for v and w as the line of symmetry y = z may not be crossed. 

For large Reynolds number, the asymptotic character of the main stabil- 
ity equation is set by a simple turning point and the resulting critical 
layer. Uniformly valid approximations to solutions of the main stability 
equation are derived using generalized Airy functions. Key steps to achiev- 
ing uniformity on the full semi-infinite domain include definition of an 
appropriate Langer variable and reformul ation of the equation in terms of a 
a proper large parameter which suitably incorporates the limiting behavior 
of the mean velocity components. Approximations to the balanced solution 



are found to have a phase shift across the critical layer which is an indi- 
cation of instability. 

A paper describing this research has been accepted by the Royal Society 
of London for publication in the Proceedings. It is co-authored by Or. M. 

Y. Hussaini of the Institute for Computer Applications in Science and Engin- 
eering (ICASE). Publication date will be late 1984. The paper was communi- 
cated to the Royal Society by Prof. Keith Stewartson shortly before his 
death. A copy of the manuscript is included in this report as an Appendix. 

The major portion of this research was carried out under the present 
grant. However, initial work on the corner layer stability problem was done 
at ICASE prior to the starting date of the grant, and appeared in an ICASE 
report. Final details of the work, including the final revisions of the 
manuscript, were done on an unfunded basis after the expiration date of the 
grant. 
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M. Y. Hussain! 

Insclcuce for Computer Applications In Science and Fngtneerlne 

Abstract 

TMs work examines Che staMllty of viscous, Incomnressihl e flow along a 
streamwise corner, often called the corner boundary layer problem. The semi- 
infinite boundary value problem satisfied by snail-amplitude .disturbances in 
Che "blending boundary layer" region is obtained. The mean secondary flow 
Induced by Che corner exhibits a flow reversal in this region. Uniformly 
valid "first approximations" to solutions of Che governing differencial equa- 
Clons are derived. L’niformicy at infinity is achieved by a suitable choice of 
the large parameter and use of an appropriate Langer variable. Approximations 
to solutions of balanced type have a phase shift across the critical layer 
which is associated with instabilities in the case of two-dimensional boundary 


layer profiles. 



IKTRODDCTIOM 


Viscous, Incompiresslble laminar flow along a streamwlse corner formed by 
the Intersection of two solid surfaces Is often referred to as the corner 
boundary layer problem. In the theoretical model of this problem, two flat 
plates, which are Infinite In the streamwlse and spanwlse directions, 
intersect such that the line of Intersection (the corner line) is normal to 
the plane containing their leading edges. The basic flow Is along the line of 
Intersection. Most theoretical work assumes a rectangular corner. This 
situation Is shown in Flg^ire 1. 

Solution of the exact corner flow velocity profile from the full Navier- 
Stokes equations is Impractical, and the only exact solution produced to date 
[10] requires an exponentially increasing suction through the walls. 
Consequently, most theoretical work on the corner layer problem has been 
numerical or asymptotic in nature with the main focus being on the boundary 
layer regions close to the solid walls. An excellent summary of existing work 
has been given by Zamir [11]. The qualitative features of the flow field 
which emerge are shown in Figure 2. Region I, away from the solid boundaries. 
Is a basically inviscid region. To lowest order, the flow in this region is 
Irrotatlonal flow. Region II, which lies "close” to the corner line 
y “ z “ 0, Is a viscous corner layer dominated by the contralnts of no-slip on 
the walls and symmetry on the line y z. Regions HI and IV are the 
boundary layers on the solid walls "away" from the corner. By symmetry, it is 
only necessary to consider region III (say) which Is the boundary layer on the 
plate y = 0, X > 0, z > 0. If the mean velocity field Is (u,v,w), then In 
region III, u Is the streamwlse component, v is the outflow from the 
boundary layer, and w Is the secondary flow component due to the presence of 



the corne'r and the opposite wall. There has been some disagreement over the 
nature of the secondary flow. However, recent experimental results support 
the picture derived by Rubin [8), Rubin and Grossman {9) and Ghla [2] which 
Involves a reversal of the secondary flow component for 0 <, z/x « 1 in 
region III. In this regard, the corner boundary layer problem Is similar to 
the boundary layer on a wing along a wing~body junction. For large z In 
region III (far away from the corner), the usual Blaslus profile must be 
recovered. However, the decay of the secondary flow with Increasing z 

appears relatively slow, algebraic rather than exponential. There will thus 
be a significant portion of region III where corner effects are Important. It 
will be convenient to denote this region where the corner layer blends Into 
the Blaslus'type layer as the "blending boundary layer." 

Theoretical work on the corner layer problem has primarily considered the 
Ideal situation where the Incoming uniform stream is aligned with the corner 
line. Hence, the pressure gradient along the corner line Is zero. However, 
it Is difficult to produce experimentally a laminar corner flow with a nonzero 
angle of incidence between the Incoming uniform stream and the corner line. 
'Chls fact, as well as apparent differences between the various experimental 
studies, has recently been analyzed by Zamir fll’. Based on his experimental 
results, v;hlch indicate a loss of profile similarity, he concludes that 
differences are the result of instability of the mean corner flow velocity 
profile. By extrapolating his data, Zamir predicts flow instability at zero 
angle of incidence when the Reynolds number is order 10^, where 

U Is the local free stream velocity, x* is the distance from the leading 
edge, and v is the kinematic viscosity. 
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The present work is the first part of an asymptotic and numerical 
investigation of the stability of the corner boundary layer problem to small 
amplitude disturbances. Three-dimensional velocity perturbations are Imposed 
on the mean flow in the portion of region III where z/k I9 small. The 
differential equations and boundary conditions governing the amplitude of 
disturbance qualities are derived, and, for high Reynolds number, the 
resulting eigenvalue problem Is found to be of singular perturbation-type. 
From an asymptotic point of view, the character of this singular perturbation 
problem is determined by the simple turning point of one of the differential 
equations and the resulting critical layer in the blending boundary layer 
region. Both heuristic and uniformly valid asymptotic approximations to 
solutions of this disturbance equation are derived. One of the two Invlscid- 
type solutions Is found to behave logarithmically avjay from the turning point 
leading to a phase shift in the general solution across the critical layer. 
Phase shifts of this type are characteristic of the critical layer in the case 
of two-dimensional boundary layer profiles, such as the Blasius boundary 
layer. 


2. TlfE DISTOEBMCE EQDATIONS 

It is convenient to work In terms of dimensionless variables to derive Che 
boundary value problem satisfied by small amplitude disturbances to the mean 
corner layer flow. Accordingly, velocities will be scaled relative to the 
local free stream velocity U. A characteristic length L may be based on 
either x* or the boundary layer thickness 6*. The choice adopted here is 
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In which case the usual Reynolds number becomes 



R 1/2. 


( 2 . 1 ) 


( 2 . 2 ) 


We wish to Impose perturbations on the mean velocity and pressure fields 
in the portion of region III closest to the corner where 0 < z/x « 1. This 
subregion will be referred to as Ilia. To lowest order, the principal 
difference between the mean velocity profile In region Ilia and the usual flat 
plate Blaslus profile Is the secondary velocity component w which changes 
direction within the boundary layer. As noted, for z/x large, w decays 
algebraically so for fixed x the Blaslus profile Is recovered as z + 

For z/x small, however, the algebraically decaying portion of w enters at 
higher order. 

Asymptotic approximations to the mean velocity profile In region Ilia have 
been derived by Rubin [8], The expansions have the forms 

u = Uq + R ^ u^ + 0(R 2). 

v“R^Vj^+0(r2)^ 

and 

w “ R~^ Wj^ + 0 (r 2^ 5 — 

X 

where ^g(s) and v^(s) are the usual Blaslus quantities 


(2.3) 
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Uq(3) - r(s), v^(3) - 2"^^^ [sf'(s) - f(s)). (2.4) 

In the current scaling, the slmilarty variable s =■ 2^^^ y, 


f'" + ff" ^ 0, 

f(0) « f'(0) =• 0, f'(s) 1 as s + 


(2.5) 



H5" + (f H5)'= 1 

H'(0) - 0, H"(0) - H' > 1 as s > 


(2.7) 


The behavior H' Is shown In Figure 3. In particular, H'(s) Is negative for 
0 < s < 2.15 and positive above this level. 

Consider now the order of magnitude of z In region Ilia. As In Ov:r 
scaling X = x*/L = R, If z =* 0 (rP) then z/x small requires p < 1. Also, 
we must have y < z to avoid crossing the line of symmetry, and as y may 
get large this requires p > 0. Rubin [B] has shown that In region Ilia, the 
function Uj^ is z/x times an order one function of s. Consequently, 

R ^ =■ 0(r . Similarly, the second term in w is 0(p ^,R . In 

particular, If p is the range 0 < p < 1/3, then the mean velocity profile 


In region Ilia Is 










jwjew! 




(u, Tr, w) 
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r”^ Wj) + 0(r“‘’) 


( 2 . 8 ) 


with q > 4/3. Thus, to the order required by the present theory for uniform 
"first approximations", the mean velocity in region Ilia Involves only y and 
the stability problem for small amplitude disturbances can be treated using a 
quasi-parallel flow approximation. 

Dropping the subscripts on mean flow quantities, consider now total (mean 
plus perturbation) velocity and pressure fields in region Ilia of the form 


u « u(y) + u'(x,y,z,t) 


V ” R ^ v(y) + R ^ v'(x,y,z,t) 


w •» R ^ w(y) + R ^ w'(x,y,z,t) 


(2.9) 


and 


p(y) + R ^ p'(x,y,z,t). / 


Relations (2.9) involve different velocity scales for the perturbations u', 
v', and w'. This Is similar to the treatment of Cortler vortices by Floryan 
and Sarlc [1] and reflects the fact that weak streamwlse vorticlty 
perturbations may produce large defects in the mean streamwlse flow. The 
linear theory for the Gortler problem has recently been placed on a firm 
footing by Kail [3], (4). In particular, for the Gortler problem it was shown 
that there is a wavenumber regime for which parallel flow gives a self- 
consistent approximation to the stability equations. 
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To obtain the differential equations for the small amplitude disturbances, 
the ejqjresslons (2.9) for the velocity and pressure may be substituted Into 
the Navler-Stokes equations and the resulting equations are consistently 
linearized about the mean flow. This gives 


3u* . — 3u^ 4 . i: ~ , 5 iil + Jil ^ 4 - i ~ liL 
3l “ 3x R '' Ty~ J~ dy R '' 3z 


1 3p' 1 _2 . 

-I 5-r- + R ^ • 


( 2 . 10 ) 


3v'' 

3t 


+ u 


3v' 1 — 3v' 1 , dv . 1 — 3v 

J. ^ V T- + -X w 


5F"^R''57~‘^R'''d7'^R 


dz 


i. 3 p: + i 72 V' 

R 3y R 


( 2 . 11 ) 


3w' _ 9w' I - 3w' 

- — + u _ — + _ V -I— - 

3t 3x R 3y 


^ 1 , dw ^ 1 - 3w' 

R dy R 3z 


R 3z R 


( 2 . 12 ) 


and 


3u' , 1 3v' . 1 3w* « 

3lT'"R3r R9^“°’ 


(2.13) 


2 32 32 3^ du 

where V " — ^ + — «■ + — ^ Is the usua' Laplaclan and -p- 

^ ^ 3z 2 dy 


~-l/2 du 
2 di • 


This is a set of partial differential equations. However, In the quasi- 
parallel flow approximation, dependence on x, z, and t can be separated out 
by assuming normal mode solutions of the form 

(u', v', w', p') » (ii(y), v(y), w(y), p(y )} e^ (2.1A) 

Equations (2.10) - (2.13) now become the set of four coupled ordinary 

differencial equations 


h “ 


la 

R 


p + u' V , 


(2.15) 
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L- V - v' V - D p , 


L. w " + m' V , 


au-.^(l Dv- 6 w) "O 


where D « d/dy, L 2 Is the differential operator 


-a'~ 6 )-vD~ iR(au w - w) , 


and a prime on a mean flow quantity denotes a y-derivative , e.g., u' - du/dy. 

Equations (2.15) - (2.19) are usually either left in their coupled form or 
rewritten in terms of vortlclty components. For the present purposes, 
however, we wish to retain primitive variables. After some manipulation, 
u. w, and p can be eliminated from (2.15) - (2.18) to give a fourth-order 
uncoupled equation for v(y) of the form 


[dL^ D - (a +8 )l. 2 ]v + iP(au' + w')Dv 


+ 1R[ou" + I w" - i (a^ + B^)v']v = 0. 


( 2 . 20 ) 


The perturbation velocity w satisfies the inhomogeneous second-order (in w) 
differential equation involving v of the form 


(o^ + 0 ^)L 2 w - i 6 L 2 Pv - aS( 0 u' - ^ w')v , 


( 2 . 21 ) 




lUilMUMVgnn.Hf.l 
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and from (2.18) 
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u 


_1 

oR 


(iDv - 6w) . 


( 2 . 22 ) 


The scenario for determining the perturbation fields Is now clear. Equation 

(2.20) , with appropriate boundary conditions on v, Is a fourth-order eigen- 
value problem. For spatial stability, the frequency w Is real and fixed and 
the eigenvalue Is the wave number vector (a,B). For temporal stability, a 
and 8 are real and fixed and the eigenvalues Is the wave speed 

c “ u)/(a^ + 8^)^^^. Having solved (2.20) and obtained Wx)» w(y) is now 
found by solving the second-order inhomogeneous problem (2.21). Then u(y) 
can be obtained algebraically from (2.22), and p(y) is obtained 

algebraically from either (2.15) or (2.17). 

The basic equation to be solved in the current stability problem is thus 

(2.20) . Some simplification of this equation Is possible as 

DL 2 » L,R - [v' r + lR(au' + | w')] , (2.23) 

In particular, (2.20) can be written in the form 

{(0“ - - 8^) - V D - (18w + v')} (D^ - -8“)v + 18w" v 

- laR{(u - ~) (D^ - - 6^)v - u" v} . (2.24) 

Equations (2.24) and (2.21) must be complimented by four boundary 
conditions on v(y) and two boundary conditions on w(y). Mo-sllp conditions 
on the solid wall y = 0 require 
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u(0) “■ v(0) “ w(0) ■ 0. 


(2.25a) 


The continuity equation (2.18) now gives the additional condition 


Dv(0) ■ 0. 


(2.25b) 


Equations (2.25a,b) provide three of the required six boundary conditions for 
v(y) and w(y). 

For modes In the discrete specturra, the perturbation auantltles u', v', 
and w' must tend to zero In the free stream, l.e., as y and z tend to 
Infinity. However, the line of symmetry y H z may not be crossed. Care 
must thus be taken In translating this 'outer' behavior of u', v', and w' 
Into conditions on u, v, and w as y tends to Infinity. Fortunately, the 
restrictions that y < z to remain below the line of symmetry and z = 0 (rP) 
with 0 < p < 1/3 for the quasl-parallel flow approximation imply that z 
will tend to Infinity automatically as y tends to Infinity. Appropriate 

'outer' boundary conditions for (2.24) and (2.21) are thus 

V, Dv, w 0 as y * +“ . (2.26) 

Consider now equation (2.24), which Is the basic equation In the stability 
problem. This equation has a simple turning point at y ■” y^, where 

au(y^) - w vanishes. As the wave speed c ” w/(a + 8 ) , this turning 

point is subtly different from the usual turning point in the Blaslus 

problem. Nevertheless, the critical layer associated with y^, sets the 
asymptotic character of equation (2.24) for large R. Further, as 



u' “ u'(y^) > 0, the reduced equation 
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R« V “ (u - OXP^ - - 3^)v - u" v_“ 0 (2.27) 


obtained by formally letting R » in (2.24) with 



(l + -2) 


,2 1/2 


(2.28) 


will have a regular singular point at y^ whereas y^. Is a regular point of 
the full fourth-order equation for v(y). Dealing with this spurious 
singularity provides one of the main challenges In deriving uniform approxima- 
tions to solutions of (2.24). In this regard, (2.24) is similar to the usual 
Orr-Sommerfeld equation for the stability of two-dimensional boundary layer 
profiles. 

In treating equation (2.24), It is convenient to consider y to be a 
complex variable. Accordingly, the aim of this study is to seek "first 
approximations" to solutions of (2.24) which are uniformly valid in sectors of 
the complex plane containing the non-negative real axis. 


3. THE PROPER LARGE PARAHETER AJTD HEURISTIC APPROXIMATIONS 

For large Reynolds number, the most obvious large parameter in equation 
(2.24) is iaR. However, if approximations are required to remain valid as 

y + “, laR is not an appropriate expansion parameter. This may be most 
easily seen in connection with the WK3 approximations to viscous-type 
solutions of (2.24). It has been appreciated for some time that the W.B 
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approximations are only valid away from the turning point In certain 

restricted sectors of the complex y-plane. However, examination of the 
*• 1/2 ^ 

O(aR) term In the Poincare series portion, of the usual UKB 

approximations shows that, In the present context, these approximations also 

1/2 

lose validity for y > 0(oR) . This breakdown at Infinity may he traced 

to the use of laR as the large parameter and may be eliminated through 
choice of the proper large parameter. 

As y + +“>, equation (2.24) becomes a constant coefficient fourth-order 
equation which can be written in the factored form 

{O^ - - 6^H(D - - Y^}v » 0 (3.1) 

where ^ 

2 22''®- 

Y - iaR(l - 0) + a + 0 +-5 + 10 w_^ (3.2) 


quality 


with 


re Y > 0 


(3.3) 


and 

V = w =» 6-. (3.4) 

Thus, as y ■* •+«>, Invlsc.id type solutions of (2,24) will have the asymptotic 
behavior 

constant x exp{i (a'^ + ^2yl/2 ^ (3.5) 


and viscous type solutions will have the asymptotic behavior 
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constant x exp{± yy + — j (y - y^)} . 


(3.6) 


Equation (3.6) suggests that the proper large parameter for the corner layer 
2 

problem Is y , not loR. Accordingly, as the starting point of the present 
asymptotic theory, equation (2.24) will be rewritten In the form 


((D^ - 


g2) _ yjj _ ^ V')}(D^ - 


6^)v 


+ 16 w" V 



.tu -• 0 
i^)( | ' " i ' "Q (® 


2 -2 * u" 

a - 0 )v - v} 


(3.7) 


with 


•b 0 




+8 10 ), 


(3.8) 


and 0Q as In (2.6). 

The so-called heuristic approximations to solutions of equation (3.7) are 
not uniformly valid in full neighborhoods of the turning point y^, and have a 
number of other limitations. However, they will be considered here as they 
contain essential elements which will be needed later to construct the fully 
uniform approximations. Let these approximations be denoted by <t>(y). 
Consider first the modified 17KB approximations for equation (3.7) based on the 
largo parameter y . Application of the usual technique to (3.7) gives 

the approximations 


<t>3(y) 


Co(Y) 


I U-0 


,-5/4 -Yho(y) + fii(y) 


e 


{l - H^(y) + 0 (y"^}, 


(3.9) 




-u- 


and 
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-s/4 "^h-Cy) +h.(y) , _, , 

^4<y) - Cq(y) |u-e| e ^ {l +Y H^(y),+ 0 (y '")}. (3.10) 


where Cq(y) Is a constant which may depend on Yi 


y„ 


(3.11) 


and 


hj(y) 


l! 


V dy. 


(3.12) 


In many applications, an explicit expression for the term Hj^(y) In (3.9) 
and (3.10) Is not required. In the present context however, Hj^(y) Is an 
essential component In the uniform approximations and Is also needed to check 
that remain valid as y •*• -?«. In the usual WKB 
procedure, H^(y) Is determined algebraically after considerable manipulation. 
An alternate procedure Is to take (j>^ (say) In the form 


_ -5/4 

® H(y,Y) 


(3.13) 


with 


H(y,Y) = H(y) + y"^ H^(y) + 0(y 


(3.U) 


Substitution of (3.13) and (3.14) Into (3.7) now gives explicit expressions 
for H'(y) and H'(y). In particular, H'(y) =• 0 so H(y) H 1. The 

equation for H'(y) is 
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«;<>-) - - [| h; '-51 ’ 

—2 

+ y [i65 + "j' + ij- + <j^ + ■ ( 7 1 'i; '’1 + >’ 6 "l''o'^ 

- I "S'" "s' ■ 


This equation immediately shows H'(y) ■*•0 as y > » so Hj^(y) ~ constant 
for large y and uniformity at infinity is preserved. Solving for Hj^(y) 
gives 


H^(y) 


an 


17 u' 

"F TT 


/u-0 1 - 3/2 
(■Tzg-J 



('U-01-1/2 K fU-Q\l/2 

■ I ^ W 


+ i^Cy), 


(3.16) 


vrlth 


\iy) 


f 1 c u'i2 cu-01-5/2 Ifii" 
'IT '■T^^ “ IL-fre 


i' vi(\}-Q\-3/2 


-2 


.Itjo— . 2 . o2ifu-0i-l/2i , 

+ -j[i6w + v +— j + a +0 Jl yrg'J I • 


(3.17) 


Following the usual practice, the lower limit of integration in (3.17) has 
been taken to be y^. Also, Hj^(y) has been normalized so that its expai.sion 
relative to y^ contains no constant term. 

The approximations (3.9) and (3.10) are valid for 0 < ly-y^l < “ in 
restricted sectors of the complex plane. For example, in the complete sense 
of Watson [6], the approximation must be restricted to the range 


/ 
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To appYoximate solutions of (3.7) of inviscld type, ^(y) may be expanded 
2 

:ln inverse powers of y » e.g., 


4' “ ^^^(y) + y ^ it> ^^^(y) + 0(y“^). 


(3.18) 


First outer approximations to two solutions <j)j and if >2 of (3.7) then 
satisfy the reduced equation 


- 0 . 


(3.19) 


To express the solution of (3.19) which remains bounded as y + +<», It Is 
convenient to define the Inviscld combination 




(3.20) 


The constant Is determined by the condition that ^ behaves 

asymptotically like exp(-(a^ + 6^) '^^y} as y + ®. For bounded ly-y^l* 

expansions for obtained by the Frobenlus method. 

They have the forms 


'i)/°^(y) = (y - y^) P^Cy) 


(3.21) 


- ? 2 (y) +^^^<0)(y) In (y-y^), 


(3.22) 




where Pj_(y) and P 2 (y) are power series in (y-y^.), ^ifyc) “ 
and subscript ”c" above denotes evaluation at y^., e.g., 

u” »• u"(y^) and u'" »■ The power series . PjCy) and ? 2 Cy) will 

not be given explicitly as they are similar to the series for the solutions of 
the Rayleigh equation In the usual Orr-Sommerfeld problem. However, It Is 

again worth commenting that u', u", and ®re different here than In the 

Orr-Sommerfeld case as the definition of the turning point y^ Involves 
0 rather than the wave speed c. 

The approximation (3.21) for the "regular invlscid" solution ijij^Cy) is 

valid In the entire complex plane. However, because of its logarithmic 

behavior, the approximation (3.22) for the "singular Invlscid solution" <j> 2 (y) 
Is only valid away from the turning point In the restricted sector 
-TT < ph(Yhp(y)) < -rr/3. It should be noted that the logarithm In iji^^^y) 
Induces a phase shift In ‘&^^^(y) across the critical layer. This behavior 
Is associated with Instability In problems Involving two-dimensional boundary 
layer profiles. 


A PRELmiNARY T8AMSF0RMATI0U 

As a first step in deriving uniformily valid approximations to solutions 
of the disturbance equation (3.7), It is convenient to explicitly bring out 
the turning point nature of the equation th»-ough a preliminary transformation 

v(y) = x(n), where n = dy}^'^^. 


(4.1) 
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The Langer variable defined in (4.1) is appropriate for the present case of an 
unbounded domain, and differs from the usual bounded domain variable which is 

the present case 


normalized by u' rather than 1-0. In 


n(y) = n'(y-y^) n" (y-y^)^ + .. 


(4.2) 


for bounded ly-y^l where 


"“c " ^1^^^^^ n" - ^ 




(4.3) 


Also, 


h(y) ~ (y y}^^^ as y -t- -f«> , 


(4.4) 


-2/3 

^ “ Y , equation (3.7) now becomes 


'5 (x + fo X - (n + 6^ f Jx" - (gj + 6^ g^)x' - 


X * 0, 


(4.5) 


where 
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” 6u - V n 


- --1 


"i 


f^Cn) » ~(4u' + iiu^ + <n) + 3u V n'“^ + [lo} + 26 ^ .+ ipw + v']n'“^, 


8j(n) “ n M, 


A 

8j^(n) = -(m" + 7uu' + 6g^ - K n li) + (2rt^ + 26^ + ipw + v')u n'"^ 


(4.6) 


+ v ( m ' + 2u^ - (a^ + 

8o(n) =• -{5u + 6nu^ + 2nu' + (ct^ + 6^)nn'"^}, 

A A 2 

go(O) = -o'" 1(0 + 8^)(a^ + g2 ^ igli + + ig-..} _ ^ 


and 


p(n) » n"(y)/[n'(y)]^ . 


(4.7) 


The function u(n) has the behavior 


Mr, = u(0) 


1 “c (-’'c ^-1/3 
3 z: ly-^F 

u 

c 


(r^) and ii(n) y n"^ as n + «. 


(4.8) 


Consider next the various heuristic approximations to solutions of 
equation (4.5). For outer approximations to solutions of Invlscld type, 
x(n) may be expanded In the form 


x(8) X ^^\n) + 6^ X ^^^(n) + 0(6^. 


''4.9) 


J 
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Then, x satisfies the transformed reduced equation 


«2 X - 0 where Rj » nD^ - Rj^(n)D.- jj^Cn) 


with D ” d/dn. Solutions of (4.10) may now be expressed 


“ nQ^(n) 


X2^*^^(o) ” Q2<’1) + 5 uq In n 


where Oj^(n) and Q 2 (n) are power series In n, Q^CO) = Q2(0) » 1 and 
QjCO) = 0. For numerical purposes It Is easier to compute and 


A (0), 


r ( 0 ), 


x/^\n) - n' 


r CO), 


X2""(n) « «2"''^>') +Z: In n' 


For later reference, Q 2 (n) also satisfies the Inhomogeneous equation 


«2 02(n) “ -5Up(?.n Q' + {gg(n) + 1)Q^] 


(4.14) 


The \mi approximations to viscous type solutions of equation (4.5) are 
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Xa(n) - dg(«) Inn'^r^^^ {l - F^(n)6^^^ + 0(6^} 


«nd 


(4.15) 


- / ^ J I -2|-5/4 +2/3 ^ /, ^ r ; n«3/2 ^ 

X^(n) " dg(6) Inn I e {I + F^(n)6 + 0(6 )| 


where Fj^(n) “ H,(y). For later use, these approximations will be normalized 


by choosing the constant dg(5) to be 


d„W) - ^ .- r,;5«. 


(4.16) 


In discussing solutions of equation (4.5), we wish to exploit as iRich as 
possible certain symmetries which they exhibit In the complex plane, and for 
this purpose we consider the Stokes and anti-Stokes lines associated with the 
WKB approximations (4.15). They are defined by the conditions 


Im Y n^^^ “ 0 and Re Y n^^^ “ 0, 


(4.17) 


respectively. There are thus, three distinct Stokes and anlt-Stokes lines 

which radiate out from the turning point n ” 0 with equal angular spacing.; 

of 2 it/ 3. For marginal stability, c is real, and for modes in discrete 

spectrum < << aR(l-O). In this case, ph Y " and the Stokes and anti- 

Stokes lines in the n-plane have the arrangement in Figure 4. In the y-plane, 

the behavior of the Stokes and antl-Stokos lines is more complicated. As 
2 3/2 

hgly) - ■j h , the definitions corresponding to (4.17) are 




Im Y ‘’^(y) ■" 0 and Re Y 


(4.18) 
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These are aj;aln straight lines with equal angular spacing of 2ir/3 close to 
the turning point y^. However, away from y^, the global behavior of the 
Stokes and antl-Stokes lines In the y-plane depends on both the details of the 
Blaslus velocity profile and the specific value of 9 chosen. These 
difficulties In the y-plane provide a compelling reason for use of the 

appropriate Langer variable. 

Consider now seven exact solutions of equation (4.5) which will be denoted 
by UqCh), Uj^(n), and Vj^(n), k«l,2,3. These solutions can be uniquely 
defined (to within luiltlplicatlve factors and modulo an arbitrary additive 
multiple of Uq In the case of U^) in terms of their asymptotic properties. 
Thus, we require that Uq be well-balanced, that U|^ be (purely) balanced 
In T|^, and that Vj^ be recessive In Sj^, where S|^ and T|^ are sectors of 

the n-plane shown In Figure 4. Associated with these seven exact solutions 

will be three exact connection formulae which will not be needed for the 

present purposes. 

Uniformlly valid asymptotic approximations for these seven solutions will 
Involve slowly-varying coefficients tines rapidly-varying generalized Airy 
functions of a stretched variable ([5], [6]) 

K » n/6. (4.19) 

Proper Indices p and q for the Airy functions may be Inferred by expanding 
coefficients In (4.5) about n •• 0 and deriving Inner approximations x(K) • 
First inner approximations satisfy the equation 

AD^ X ” 0, 


(4.20) 





iVPT**IM?W 




¥ ' 
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where D » d/dK and A «= D - 5 . Solutions of (4.20) are 8^(5, 2), 

and Aj^(5,2), k-1,2,3. Higher Inner approximations are solutions 
of Inhomogeneous versions of (4.20) and Involve the Airy functions Bq(£;,p), 
®lt(C,p,l), and A|^(C,p) with p- 0 , 1 . 


5. SOLUTIOHS OP DOJCEMMTr-RECESSIVS TWE 

Uniformly valid "first approximations" to the dominant-recessive solutions 
Vj^(n) of equation (4.5) involve the Airy functions A^(C,p), k»l,2,3, with 
p“0,l,2. They have the forms 

V^(n) ~ a(n)A^(C,2) + 5 ^ h(n)A^(^,l) + 6c(n)A^(C ,n) , (5.1) 

where a, b, and c are regular at the turning point. On substituting (5.1) 
Into equation (4.5), the slowly varying coefficients a(n) and c(n) are 
found to satisfy the equations 

R2(na) => 0 2) 

and 

{2nD + (nfjj “ gj^)} (a + nc) »> 0, (5.3) 

where R 2 Is given In (4.10), D » d/dn. These equations may immediately be 
integrated to give the regular solutions 

s(n) « a( 0 ) n"’' X|L*'°^(n) » a( 0 ) Oj^(n) 


(5.4) 
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a(n) + nc(n) - a(0)(^)"^^^ exp{h^(y)} 


(5.5) 


The differential equation satisfied hy b(n) is quite Involved and will not 
be given explicitly. A more Illuminating way of obtaining this coefficient is 
to require that the outer ejcpanslon of V,^(n) (say) be the WKB approximation 
X 2 (n) away from the turning point in the sector where this solution Is 

purely recessive. With normalized as In (4.15) and (4.16), this 

condition gives 


a(0) » 1 


(5.6) 


b(n) - 2n~^ c(n) + n"^^^(a(n) + nc(n))(Fj(n) - ] . (5.7) 


Thus, b(n) Is essentially a regularized form of Fj^(n) " H^(y). 


6. SOLUTIONS OF BALANCED AND WELL-BALANCED TYPE 

The uniform first approximation to the well-balanced solution Is simply 


U()(n) - x/°^n) + 0(6^, 


( 6 . 1 ) 


First approximations to the three solutions of balanced type have the forms 





MW lilJti j. ww,>i' I'l ' '! '-I w t l'»MlfWffiMWW3CS^5f!5^^ 
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U^(n) ~ G(n) + X6(A(n) B^(5,2,l) + 6^ B(n)Bj^(C ,1,1) 

+ 6C(n) Bj^a-,0,1)}, (6.2) 

where G, A, B, and C are regular at the turning point and X Is a constant 
to be determined. 

Substitution of (6.2) into (4.5) shows that A, B, and C satisfy the 
same differencial equations as a, b, and c. Consequently 

{A(n),B(n),C(n)} - A(0){a(n),b(n),c(n)} . (6.3) 

The well-balanced part G(n) of (6.2) Is found to satisfy the equation 

G - x{4n(A + nC)' -t- A - 2n(nC)' + f^(A + nC) - nc} (6.4) 

which, on simplification, becomes 

R 2 G ■= XA(0)(2nQ' + (I + go)Qj. (6.5) 

A comparison of this equation with (4.14) now shows that the solution of (6.5) 
which Is regular at n ” 0 Is of the form 

G(n) - Q2(n) + bQ X^^°^(n), (6.6) 


where bQ Is an arbitrary constant, provided 








W5®BP’’(>**iW*WI!W! 
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(6.7) 


The constant bg In (6.6) may be chosen In a variety of ways. A particularly 

convenient way Is to choose bg so that the outer expansion of U 2 (n) (say) 

to 0(6^) In the sector T 2 contains no multiple of i^^^^^y). Using the 

asymptotic expansions of the B-type Airy functions In T 2 and the usual outer 

2 

expansion operator E 2 to 0(6 ), 

E 2 U 2 (n) - Q 2 (n) + + 5Y(,[log n - log 6 + y - 1 - 2irl]} (6,8) 


where Y “ 0.577 ••• Is Euler's constant. This result can now be written In 
the form 

^2 U^(ri) «• X^^^^^y^bg - 5Ug[log 6 - log n' + 1 - Y + 2tt1]} . 

(6.9) 

Thus, E 2 U 2 (n) contains no multiple of 4>j;^^^(y) provided 

A 

bo ■ 5Uo[lo8 6 - log n' - Y + 1 + 2nl] . (6.10) 

It is worth noting that the slowly varying coefficients In both (5.1) and 
(6,2) are expressible in terms of known quantities in the heuristic 
approximations: the regular solution and the regular part of the singular 
solution of the reduced equation, and the terms axp(hj^(y)) and H|(y) in 
the WKB approximations. 




7. CONCLUDING REMASKS 

The major alms of this paper have been to obtain the boundary value 
problem associated with small amplitude disturbances in the blending boundary 
layer region, and to derive "first approximations" to solutions of the 
governing equation for v(y) which are uniformly v-lid on a semi-infinite 
interval. This requires the Introduction of a new large parameter and an 
appropriate Langer variable. The uniform "first approximations' Involve 
slowly varying, analytic coefficients times rapidly varying generalized Airy 
functions of a stretched Langer variable. Because of the Important first and 
third derivative terms in the fourth-order disturbance equation, as well as 
the presence of v and w in other terms, these analytic coefficients differ 
from their counterparts in the usual stability problem for two dimensional 
mean profiles. Approximations to solutions of balanced type display a phase 
shift across the critical layer that is an indication of instability. 

A derivation of the eigenvalue relation associated with the boundary 
conditions (2.25) and (2.26), numerical results for the solution of equation 
(2.21) for w, and stability characteristics of the basic corner flow will be 
given in a subsequent paper.' The characteristic equation for this problem is 
of the form A(0) = 0 where 

A(y) - n; W($,6Vj)(n), (7.1) 

W is the usual Wronsklan, 

«(q) » A^®^UQ(n) + U 2 (n) - «(y), 


(7.2) 
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and Is again deCermlned by the condition that iSCn) remains bounded as 

n ♦ +«. The need to use uniformly valid approximations for Uq, U 2 , and 

In analyzing (7,1) may be easily seen. If, following the usual heuristic 
procedure, Frobenlus-type approximations are used for and either first 

WKB approximations or first Inner approximations are used for V^, the 

resulting crude approximation to (7.1) Is Independent of both the mean outflow 
velocity V and the mean secondary flow component w. 
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